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Abstract 

New theoretical formulae of the sound velocity and the B/A nonlin- 
earity parameter for some fluids are presented in the paper. Semi-ideal 
and van der Waals models of gas are considered and the parameters are 
compared to experiment data. The liquid water model for equation of 
state, given by Jeffery and Austin analytical equation, is considered also 
and modified on the basis of acoustic data . 



1 



NOTATIONS: 



x space coordinate [m]; 

t time [s]; 

p density [kg/m 3 ]; 

p pressure [N/m 2 ]; 

v velocity [m/s]; 

r], £ viscosity parameters; 

X heat conductivity parameter; 

a;*,£*,/9*,p*, dimcnsionless variables; 

T absolute temperature [if]; 

u internal energy per unit mass [J/%]; 

Po , Po i w o , eo , ?o unperturbed values; 

p 7 p,v,e,T perturbations; 

/3 characteristic scale of disturbance; 

a coefficient responsible for amplitude of acoustic wave; 

D1..D5 dimcnsionless coefficients in evolution equations; 

E1..E5 coefficients in caloric equation of state; 

c linear sound velocity [m/s]; 

B/A, C I A acoustic parameters of nonlinearity; 

s entropy [J • feg -1 • K -1 ]; 

c v ( p \ heat capacity under constant pressure (volume) [J • kg^ 1 ■ K~ 

R the universal gas constant [J • mol^ 1 ■ K -1 ]; 

Ri individual gas constant [J • — 1) ■ kg^ 1 )] 

fi molar mass [kg/mol]; 

f osc number of oscillation degrees of freedom of a gas molecule; 

0i characteristic temperature of oscillation [K] ; 

7 adiabatic gas constant (c p /c v ); 

a, b van der Waals constants; 
X,ai,v B ,T B , 
a, b , h, b 2 , 

Aq, ^1,^2, ^3 constans in Jeffery- Austin equation for water; 



1 Introduction 



The experimental researches of some physical properties of different fluids, such 
the sound velocity and the acoustic nonlinearity parameters, are well known 
and advanced today ^lElElEl- The theoretical basis on these problems, still 
have a lot of aspects to be studied properly. 

One of the most interesting thing is the connecting of statistical thermo- 
dynamics and acoustic studies for gas or liquid, in order to disclose the micro- 
properties of a nonlinear propagation medium in a direct link to the macroscopic 
one. For example, knowing the structure on the molecular level for a model 
medium and comparing its acoustic properties to real one, we could conclude 
about a molecular structure of the fluid. In order to obtain the information we 
do not have to solve the system of basic equations. Obviously, we need both 
thermic and caloric equations of state, taking into account the thermodynamical 
relations between them 5 . 

There are used two different representations of the equations of state in our 
paper: the Taylor series for thermodynamic variables in a vicinity of a mechan- 
ical equilibrium point 00?] and some analytical formulas. We start from 
well-known ideal/semi-ideal and van der Waals gases, comparing and discussing 
the results of the sound velocity and B/A evaluations in both approaches and 
experiment. Hereby we consider an analytical (thermic) equation of state for a 
liquid water accordingly to a formula given recently by Jeffery and Austin [Hj. 
Application of this equation to find sound velocity c and B/A is realized taking 
into account the equation for free energy. We presents a general form of the 
formula for the sound velocity and the nonlinear parameter B /A, (compare with 
the popular formula p^) with and without using the mentioned Taylor series. 
(The theme of higher order parameters of fluid is presented in , look also in 

HH. ) 

The other interesting question is using the sound velocity, and nonlinear 
parameters, to test some new equations of state. We mean that an experimental 
value of the sound velocity (and B/A, C/A...) can be compare with a theoretical 
one, so then we can except or accept, a new model of medium. More, we 
have a mechanism to make some corrections, first of all we mean the adjusting 
parameters choice whilst covering more vide field of applications . 

The following section includes formulating the physical problem on the math- 
ematical level, similar like in [?, 7 , with using projecting technique. Widely, 
this theme was raised earlier in jjjj. Next sections contain general approach to 
the fluids parameters and its adaptation to some individual theoretical models. 

2 Formulating of mathematical problem 

The considered physical problem is the fluid medium (gas or liquid) being under 
acting the acoustic wave. A basic system of the hydrodynamic laws of conserva- 
tion of momentum, energy and mass in one-dimensional How is given by known 



equations: 



dv ^ dv ^ 1 dp /4 _|_ A 9 2 v q 
9i 9a; p 9a; \ 3 / 9a; 2 

du du dv (A \ ,9u x2 9 2 T „ m 
'ft + PV 8~ X +P d- x -{^ + C ) fe> ~ *8* = °' (1) 
9p 9(/w) = Q 
9t 9a; 

and their simplified forms for nonviscous and non-heat-conducting fluids: 

dv dv 1 dp 
dt dx p dx 

du du dv n /„.. 

p m +pv d-x +p d-x=°> (2) 

dp d(pv) 



dt dx 



0. 



In order to complete the physical problem we use the caloric and the thermic 
equations of state. The general forms of these thermodynamic equations, ob- 
tained by the Taylor series of two variables (p, p) are: 

p u = Eip H p H p H 2— p H pp + ... (3) 

Po Po P P 

f=^p + ^p + .... (4) 

P0Cv PoCv 

Obviously, we assume the quantities u, T 7 p, p, v have to be treated as z — z D + z 
(index "zero" means an equilibrium value of z and "prime" means an addition, 
caused the sound wave), so the above formulae was written for the additions 
only. Finally, using dimensionless variables: 

v = acv*, p = ac 2 p p*, p = ap p*, x = fix*, t = t*0/c, 

where c is the linear sound velocity, given by 



[3 means the characteristic scale of disturbance along x and a is the coefficient 
responsible to the amplitude of the acoustic wave, we can formulate problem as 
the one matrix equation: 

JU + £* = * + * + 0(a 3 ), *=| /) ). ((i) 




This is the nonlinear operator equation of time evolution, where: 



L = 



dx 



dx 



d_ 

V dx 







* = a 



I dv ,dp \ 

V dx ^ dx 

dp dv . , , _ . 
9p ,dv 



(7) 



9a; 



9a; 



The asterisks for the variables were omitted for simplicity, and D\ , D 2 denote 
dimcnsionless coefficients, which are algebraic functions of E1..E5 (see 

^ = ±(-1 + 2^* + *), 

The second-order nonlinearity column will contribute to the B j A nonlinearity 
parameter. The constants A, B relate to coefficients Dj i Ej in the following 
way: 

A = [(1 - E 2 )/E 1 ]p , B = —{D\ +D 2 + 1)[(1 - E 2 )/E 1 } Po , 



I = -D 1 - D 2 - 1. (8) 

Now, let us return to the evolution equation © . The new application of method 
of acting projectors was presented also in That is the simple way of sep- 
arating the leftwards, rightwards and stationary modes of sound, which are 
responsible for the wave propagation effect in 'left' and 'right' directions, and 
for other effects (such as "streaming", seejlj). The separating of mode is done 
on the evolution equation level. Acting by one of mentioned projectors, i.e. 
unitary, orthogonal operators Pi, P 2 : or P3 



, ^3 = 



on the evolution equation gives us a new form of wave equation: 

j^i* + Pi* + P* = 0, * = 1,2,3, (9) 

and introduces the sound velocity c. The simplified version of the equation, 
neglecting heat conductivity and viscosity of medium, has form: 

i.m 





i = l,2,3; m = l,2,3; u=l,2,3; c\ = 1; C2 = — 1; C3 = 0; 



where T denotes some coefficients matrix, and it is built with algebraic sums of 
D\ and D 2 , for 11= 1: 



Y 1 

i = 1 
i = 2 
i = 3 



m = 1 m = 2 r7i = 3 

-£>i - -D 2 + 1 Dt + D^-l 

- D 2 - 3 Z?i + L> 2 - 1 

-L>2 - 1 D 2 - 1 



We have to remember that the below form of equation is written for dimension- 
less variables. 



3 Velocity of sound and B/A parameter in 
medium 

In order to apply our considerations for some different models of fluids, not only 
ideal gas, we use the general expression for the sound velocity. Here, the sound 
velocity c takes the part of coefficient in the wave equation. 

<?={¥) w 

\°PJ S,p=p ,p=p 

Index 'iS" means obviously that entropy is constant, however in practice it is 
enough to assume the reversible adiabatic process [1]. 

The sound velocity dependence on temperature is written well in the exper- 
imental acoustic papers. So, some experimental figures show linear dependence 
c on T for majority of liquids. However, it must be notice, that the sound veloc- 
ity grows due temperature to 74°C and next becomes smaller for liquid water 
case. These special forms of curves c(T)-dependence are known in literature as 
Willard curves. The peculiarities of water result from long-range order, strong 
polarity and strong association of water molecules. 

The previous papers providing experimental data of B /A such ^1 E] and 
later papers, show that the ratio B/A generally increases slowly with tem- 
perature, although there are some exceptions. The contribution to B/A from 
temperature changes is smaller than one due pressure changes. 



3.1 General formulae for nonlinearity parameters 

To find a formula for the sound velocity we need two equations of state: p = 
p(p,T), U — U(p,T) and the first law of thermodynamics: dU — TdS — pdV, 
where S means entropy, in form: 

du = Tds + -^rdp. 

p 2 

(s,u are variables expressed per unit mass.) The differentials dp and du we 
take as: 

dp={^) dT + {^ T dp^ M T + P2 dp 



du = (J?) dT + (fiO dp = MT + Mp 

Comparing suitably the expressions and introducing them to the first law of 
thermodynamics equation, one can obtain a formula: 

dp Pi f P a , PsPz \ , Pi^ds no , 

d- P -jA7- p * + —) + j/d- P - (12) 

Next, we can do an assumption of adiabatic process of propagating sound, so 
finally: 

„2 



where 



P: 



P: 
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( dp 
~ [dp 




Pa 


(du 
~ [dp 



(13) 



That is a new approach, without using of Taylor series for variables. 

If we use the Taylor expansion for relation between pressure and density, and 
limit ourselves to quadratic terms, we will get the expression for B/A nonlinear 
parameter in form: 



A 4 \dp 2 J POiS c\ \ dp j 



(14) 



3.2 Ideal and semi-ideal gas model 

For ideal gas we receive the coefficients in forms: E\ = E4 = 1/(7 — 1), E% = 
E$ = —1/(7 — 1), £3 = 0, D-y = —7, Z?2, so using the mentioned equation © 
we receive @]: 

hpo 

Po 

where 7 = C p /Cv- 

B/A = 7-1, C/A=( 7 -l)(7-2) 

The semi-ideal gas model accepts energy of oscillation of molecule, and omits 
energy of rotation, because it is significant for very low temperatures and light 
gases only. The model concerns polyatomic gases only, because for monoatomic 
ones we have the same formulas as before. 

(Mc v ) sld = [Mc 

v)id 

+ (Mc) osc + A(Mc) 

rot 

We use the Einstein - Planck formulae for vibrational specific heat: 

Use /n x 2 e e t/T 



(Mc) osc = (MR)<T(j;) 



by using which we get the equation for internal energy for semi - ideal gases [6] : 

where Mc, Mi?, fi, f, 9i - vibrational specific heat, universal gas constant, molar 
mass, number of degrees of freedom of a molecule and characteristic temperature 
in correspond order. According to J5j the sound velocity formula looks finally: 

c> = R.T | 1 + + £ (|) 2 (e^ - l) ^ (16) 

A formula for B/A has a more cumbersome form: 



where some new symbols mean accordingly: 

^ = -- 1 T +E{^(-c 3i -i)} 

7 — 1 *r 



C 



2 1 



tfA 2 exp(^/T ) 



T o ; (exp(^/T ) - l) 2 

exp(^/T ) 



C 3i = -2- (i%/T ) + 2{#i/T ) 
C« = (tf 4 /T 
C 5 , = -1 + 2 



(exp(^/T ) 

(exp(^/T ) - l) 2 
exp(iVT ) 



(17) 



(exp(^/T ) - 1) 

The semi-ideal gas model has provided the quite correct data for polyatomic 
gases, like C02,CH4, for any monoatomic (no oscillations) and diatomic gases 
we have not interesting difference for both models: semi-ideal and ideal one. 
Some results for a few gases are presented in Table 1 and the temperature 
dependence of c and B/A for C0 2 is shown at Fig.l. and Fig. 2. (more in 0). 



Gas 


Model of ideal gas 


Model of semi-ideal gas 


Experimental data 




c[m/s] 


c[m/s] 


c[m/s] 


He 


972.9 


972.9 


971 


co 2 


262.2 


255.0 


256.7 


CH A 


434.7 


431.3 


430 



TABLE 1. All values in the table are obtained for T = 273.15^. 
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Figure 1. 

CC>2- Comparison of temperature dependence of sound velocity 
for both theoretical models in 200-480 K range of temperature. 
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Figure la. 

CO2 and CH4. Differences of sound velocities 
Cid — c.sid for N 2 , CO and CO2 gases. 
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Figure 2. 

CO2 and CH4. Comparison of B/A for both theoretical models 
in 200-520 K range of temperature. 

3.3 Van der Waals gas model 

Van der Waals gave the famous equation of state for gas model in 1873 year. 
That is convenient to use the mentioned equation in form: 



P 



pRiT 
(1 - bp) 



ap 



and 



e = c v T — dp, 



(18) 



where Ri = R/p means individual constant for gas, a — a/p 2 , b = b/p - van der 
Waals constants and Ri, e and c v denote some values per unit mass. Calculating 
of the sound velocity c, according to the (0 gives: c 2 = 



Ri T 1 



c v a{\ — bp) c v (p + ap 2 )b c v (jp + ap 2 ) (1 — bp) 



Ri pRi 

and by the new formula . gives a formula: 

7-1 



p 2 R t 



1 



pb 



(1 -bp), 



2dp; 



a j p-p 
(19) 

(20) 



c 2 = R* T 

[l -bp) \(l-bp) 

Using JSJ) provides B/A in follow form: 
Be ' 

— = \ R 2 p + 2RiC v dbp 3 + R t c v p + dR 2 p 2 + 

+ 6c 2 v dbp 3 + 2c v R l bpp - 2c 2 d6 2 p 4 + 2c 2 v bpp+ 

+RiC v dp 2 } |(i?iP — c v dp 2 + 2c v dbp 3 + c v p + Ridp 2 )c v (—l — 6p)| 

(21) 

Table 2. contains a comparison of some results for sound velocity for 273. 15K. 
Fig. 3. presents the theoretical curve of pressure dependence for c. 



0^(1 -bp)' 1 



Gas 


Model of ideal gas 


Laplace formula 


Model of van der Waals gas 


Experimental data 




c[m/s] 


c[m/s] 


c[m/s] 


c[m/s] 


He 


970.9 


971 


970.7 


971 


H 2 


1259.9 


1261 


1259.2 


1286 



TABLE 2. All values in table are obtained for T = 273.15.fC 
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Figure 3. 

CO2. Comparison of theoretical (van der Waals model) 
and experimental values |12j of sound velocity for changing pressure. 
(C0 2 , 323.95.fY, 0.3MHz) 

The van der Waals equation of state in its standard form, give good results 
for some gases, but is not valid for many of liquids, in particular for water. 



3.4 Liquid water model 

An analytical equation of state for liquid water 

The trials of "build" an analytical equation of state for liquid water were made. 
The one of the newest is the Song-Mason-Ihm equation J2j and its modifica- 
tions, for the polar fluid, made by Jeffery and Austin [Hj: 

ap a lP 



pRT ur iff 1- Xbp 

where p is expressed in (mol/m 3 ) unit and function b(T) has form: 

b(T) = v B (o.25e 1 /( 2 - 3T / T - +a5 ) - b ie 2 - 3T ^ + 6 2 ) 

The constants used by authors of paper have the following values: A = 0.3159; 
ai = 2.U5v B ; b Q = 1.0823u B ; h = 0.02774; b 2 = 0.23578; v B = 4.1782 * 
10- 5 m 3 /mol; T B = 1408.4X; a = 0.5542Pam 6 /mol 2 . 

The equation of state for free energy proposed by Jeffery and Austin for T > 4°C 
has form: 

F = Ai(p, T) - RT^(T), (23) 



* = * 1 + * 2 ^L^ + *3^, 

where: A Q = 21A7kJ/mol, *i = 5.13, * 2 = 20.04, * 3 = 2.73, a A means 
temperature wavelength: 

A2 R 5 / 3 h 2 

2irmK 8 J 3 T' 



In order to find the acoustic wave propagation velocity according to (|13|l in 
discussed medium, we make the following calculations for adiabatic process. 
The known expression for free energy will make possible finding internal energy 
per unit mass. In statistical physics: U = F — T (fy) y , so dll = dF — dTFx — 
TdFx- Here, bottom index means partial differential We can write: 

dF = dAi - RdT^ - RT^ T dT, F T = A 1T - - RTf T , 

dF T = dA 1T - - RdT^T - i?T* TT dT, 
dAi = Ai T dT + A lp dp 
and in the same way like above we obtain the equation: 

dp (a 1p - TA 1Tp + ?jj£\ + dT {-TAitt + 2RT^ T + RTH TT ) = 

Sound velocity and B/A in liquid water 

and finally: 

(A lp - TA 1Tp + + %TA 1TT ~ ^2RT^ T - ^RT^ TT ) fa 

c 2 = 1 1 L f24) 

TA 1TT -2RT^ T - RT 2 m TT ' v ; 

The expression of the nonlinear parameter B/A, received by using (|14(l has 
more complicated form, but it is not difficult to calculate some values using a 
computer. 

Results for analytical model of liquid water 

Authors tested Jeffery-Austin analytical equation of state for liquid water 
[8], being a development of |13| . Below we present some diagrams for the 
sound velocity and the nonlinear parameter B/A. The equation seems to be 
rather sensible for small changes of constants. Some results of c and B/A for 
the Jeffery-Austin equation differ from some experimental data (for 5 — 55°C 
temperatures), what was shown on the following figures (from Fig. 4. to Fig. 7.). 
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Figure 4. 

Water. Dependence of sound velocity 
on temperature in 10 5 Pa pressure. 
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Figure 5. 

Water. Dependence of sound velocity 
on pressure, T = 303.15 K. 
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Figure 6. 

Water. Dependence of B/A parameter 
on temperature in 10 5 Pa pressure. 
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Figure 7. 

Water. Dependence of B/A parameter 
on pressure, T = 303.15 K. 

However, the most interesting and worthy to underlining is fact we have 
a mechanism that makes possible testing and even correcting the equations of 
state. One of some various possible corrections is changing the constants: A = 

0. 244 instead of 0.316, b a = -0.000026 instead of 0.000045 and /psi 2 = 22.04 
instead 20.04. The figures above shows compatibility that changed model to 
some experiments. The presented corrections have an example character only. 

4 Conclusions 

1. The new presented formula for the sound velocity is better than the earlier 
known one, used by the other authors, because of no necessary using of Taylor 
series, we have no need to limit ourselves to some first expressions in that 
expansion. Although, the values given by the both of these methods arc the 
same in the semi-ideal and van der Waals gas model cases, probably we can 
expect interesting results of comparison for more complicated models of fluids. 
Undoubtedly, this is the small step forward in the theory. 

2. Connecting thermodynamic physics and acoustics seems to be an in- 
teresting source of information about considered medium. We make a sensitive 
mechanism to test and correct theoretical models of various fluids, using of some 



experimental data for c and B/A. In future, incorporating the links between 
statistical physics and thermodynamics, it could be possible concluding about 
intcrmolecular potentials and molecular structure of medium from acoustic re- 
searches of fluids. 
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